Time-dependent self-trapping of Bose-Einstein Condensates in a double-well potential 
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Based on the mean- field approximation and the phase space analysis, we discuss the dynamics 
of Bose-Einstein condensates in a double-well potential. By applying a periodic modulation to the 
coupling between the condensates, we find the condensates can be trapped in the time-dependent 
eigenstates of the effective Hamiltonian, we refer to this effect as time-dependent self-trapping of 
BECs. A comparison of this self-trapping with the adiabatic evolution is made, finding that the 
adiabatic evolution beyond the traditional(linear) adiabatic condition can be achieved in BECs by 
manipulating the nonlinearity and the ratio of the level bias to the coupling constant. The fixed 
points for the system are calculated and discussed. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in a double- well po- 
tential have attracted much attention in the past decades, 
it provides a useful tool to study fundamental problems 
in quantum physics at the macroscopic scale and opens 
the possibility of practical applications, such as high- 
precision measurements, interferometry and thermome- 
try 11] . Different from Josephsoniunctions realized in 
superconductors or superfiuids [2|, y] , the interatomic in- 
teractions in BECs play an important role in its dynam- 
ics, leading to many rich and interesting nonlinear effects. 
Self-trapping [J-u| is one of these interesting phenomena 
from which much attention has been received in recent 
years. When the initial population imbalance between 
the two wells and the nonlinearity in BECs is above a crit- 
ical value, the amplitude of the Josephson oscillations are 
extremely compressed and most of the atoms are trapped 
in one well as theoretically predicted and experimentally 
observed Q. 

The self-trapping is originally defined and studied for 
a time- independent system (Hamiltonian). Consider a 
system governed by a time-dependent Hamiltonian, the 
following questions naturally arise: How can we make the 
self-trapping happen at a time-dependent state? What 
are the fixed points in this situation? How do these fixed 
points behave? We will answer these questions in this pa- 
per. In fact, adiabatic dynamics in BECs has been inves- 
tigated in different regimes, such as Landau-Zener tran- 
sition Rosen-Zener process [l7], EH- Meanwhile, 
several proposals have been proposed with the stimulated 
Raman adiabatic passage technique(STIRAP) [19] , to co- 
herently manipulate BECs in double-well or triple-well 
potentials by adiabatically following a spatial dark state 
of the system [20l - f2^ | . finding that nonlinearity plays a 
neg ative role against adiabatic evolution in some regimes 
[24j. Different from these works, we discuss here the pos- 
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sibility of manipulating the BECs to adiabatically follow 
the instantaneous eigenstate of the system with the help 
of atom-atom couplings. We find that nonlinearity can 
help the system to track the instantaneous eigenstates 
of the system, leading to a self-trapping. By manipulat- 
ing the nonlinearity and the ratio of the level bias to the 
coupling constant, adiabatic evolution beyond the tradi- 
tional adiabatic condition can be achieved in BECs in 
both symmetric and asymmetric double-well potentials. 

The paper is organized as follows. In Sec. UH we in- 
troduce the model and transform the equations to the 
time-dependent representation, choosing the instanta- 
neous eigenstates as the basis. In Sec. IIIII and Sec. IIV1 
we investigate the self-trapping in BECs in a symmet- 
ric and an asymmetric double- well potential respectively, 
and discuss the effect of nonlinearity on the dynamics. 
Finally, we conclude our results in Sec. [V] 

II. MODEL 

We start with the standard nonlinear two-level model, 
which describes a BEC in a double-well potential [1, [H, 
H3| . In a time-independent basis {\L),\R)} : the model 
can be written as [J, Q 

. d fa\ 1 (r + ms v \ / a \ 
l ~dl\ b ) 2\ v -r-ms)\b) 

where s = \b\ 2 — \a\ 2 stands for the population imbalance 
between the two wells denoted respectively by \L) and 
\R). r, v and m characterize the level bias, the tunneling 
between the two wells and the nonlinearity in BECs, re- 
spectively. Here and hereafter, we rescale to, r, v in units 
of w, and t in units of 1/cj, h = 1 has been set, hence all 
parameters in this paper are of dimensionless. In previ- 
ous study, the tunneling v is a real parameter, here we 
discuss a complex tunneling dependent on time through 
V = ve lujt , which can be viewed as a nonlinear two-level 
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system driven by a time-dependent magnetic field [25| . 
B = (B Xl By, B z ). There are many ways to realize ex- 
perimentally such a Hamiltonian. For instance, one may 
simulate this Hamiltonian by a BEC in a double-well po- 
tential with the height of the potential barrier modulated 
and the phases of the two modes mismatched. Define 
cos/3 = , T 5 , a = 2B — Vr 2 + v 2 , d> = Lot, the effec- 

tive Hamiltonian H (replace v by V = ve lujt and V* 
respectively in Eq.(fl])) can be rewritten as, 

Hq = Bq (sin /3 cos (f>a x + sin /3 sin (\>a y 

+ cos (3a z ) + ^ms(T z , (2) 

where Bo acts as the amplitude of the magnetic field, 
o x , <j y and a z are the Pauli matrices for the two- level 
system. /3 behaves as the ratio of the level bias to the 
coupling constant, a characterizes the amplitude of the 
driving field, and <p represents the swept angle of driving 
field. 

It is easy to write the instantaneous eigenstates of 
the system, \E+) = (cosf, sinfe'"*) 7 and \E_) = 

(— sin je~ tult , cos in the basis \L) and \R), with the 
corresponding eigenenergy E± = ±a/2. After a repre- 
sentation transformation by replacing the original time- 
independent biases \L) and \R) with the instantaneous 
energy eigenstates \E + ) and \E_), we obtain a new set 
of equation 

.d_ fc\ _ (E+ + M x s + E ++ M 2 s + E+- \ ( c\ 
% dt\d) \ M 2 *s + E_+ E_ - M x s + E— ) \d ) 

(3) 

with c and d representing the amplitude of the BECs in 
\E + ) and \E_), respectively, i.e., \ip) = c\E+) + d\E_). 
The following notations have been used, 

s = cos/3si + sin /3s2 

= cos (3(\d\ 2 - |c| 2 )+sin/3(cd*e # + c*de-^), 

Mi = — cos/3, 

M 2 = -^ S in/3e-^, 

E ++ = -E — 

= -i{E+\E+) 

= |(l- C os/3), 

E+- = E*_ + 

= -i(E + \E_) 

= |sin/3 e -^. (4) 

Here s representing the nonlinearity consists of two types, 
namely si and s 2 , which correspond to population imbal- 
ance and spatial atomic coherence (2^, respectively. Mi 
and M 2 characterize the strength of nonlinearity in the 
system. E± and E±± are actually the integrand of the 




FIG. 1: (color online) Illustration of self-trapping states on 
the Bloch sphere. The Bloch vectors for these states are 
f x — sin (27 — ft) cos </>, f y = sin (27 — /3) sin</> and f z = 
cos (27 - /3). The basis is \\L), \R)}. The red-solid lines de- 
note the instantaneous eigenstates of the system and blue- 
dashed lines represent the other self-trapping states. The pa- 
rameters chosen are r = 1, v = \/3 and /? = tv/3. 

dynamical phase and Berry phase [2^,[23- -Eh an d E (_ 

are terms that may induce tunneling between the new 
bases \E + ) and |EL ), and they are also the par ts which 
lead the adiabatic theorem to break down [28| for the 
system without nonlinearity. Observing Eq.([3]), we find 
that the nonlinearity may help to track the system on the 
instantaneous eigenstate of the Hamiltonian, meanwhile 
the change rate of the Hamiltonian can affect the self- 
trapping of the BECs on the instantaneous eigenstates, 
how do they relate to each other? 

With the mean-field approximation, the probability 
amplitudes can be written as c = |c|e ie<: and d = |d|e l6,£i . 
By defining population imbalance and relative phase as 
Z = \d\ 2 — |c| 2 and = 8d — C , respectively, we can 
study th e dy namics of the system in its (classical) phase 
space pLllOl]. However, the situation under study is little 
bit different, the BECs are driven effectively by a rotat- 
ing magnetic field with frequency uj, this motivates us to 
define pseudo fixed points by 

Z = 0, and 6 = u, (5) 

in the phase space, which corresponds to states of the 
form, = \c\\E+) + \d\e iLUt \E^) with fixed popula- 
tion difference in the language of wavefunction. We 
can track this state both in adiabatic and diabatic evo- 
lution. For simplicity, we will refer the pseudo fixed 
points as fixed points when no confusion arises. Define 
sin 7 = |c| and cos 7 = \d\, the state corresponding to the 
fixed point can be represented in the basis {\L), \R}} as 

|^) = (sin( 7 -|), cos( 7 -|)e-*) T . 

By manipulating 7 (location of the fixed points), we 
can obtain different instantaneous self-trapping states, 
which maintains the population imbalance and keeps the 
relative phase matched with the driving field, as shown 
in Fig. [T] For example, the fixed point 7 = (Z = 1) 
and 7 = \ [Z = — 1) represent the instantaneous energy 
eigenstates = \EJ) and — \E + ), respectively. 
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FIG. 2: Phase-space trajectories with different nonlinearity. 6 is in units of 7T. v = 1 and oj = 1. The nonlinearity is set 
as m — in (a), m — 2 in (b) and m — 4 in (c), respectively. For weak nonlinearity, there exist two fixed points. When 
I A/2 1 > I £7-1 — I (strong nonlinearity), the other two fixed points appear in the phase-space. The system can be trapped in the 
pole fixed point (we denote by the pole fixed point the point with Z very close to ±1, i.e., near the pole on the Bloch sphere) 
when nonlinearity is strong enough. The inset in (c) is the enlarged circles around pole fixed point at the bottom of the phase 
space. 



III. SYMMETRIC DOUBLE- WELL POTENTIAL 

In this section, we focus on the BECs trapped in a 
symmetric double-well potential (i.e., r = 0, j3 = tt/2) 
and discuss the self-trapping with different nonlinearity. 
In this case, the parameters in Eq. © reduces to 
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= &2, 


E± 


= ±v/2, 


Mi 


= 0, 
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M 2 




E++ 


= w/2, 






E+- 


= V, 
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(6) 



where only the nonlinearity s 2 from the spatial atomic 
coherence takes place. In this case, Eq. (j3j can be writ- 
ten as 



. d 

''dt 



v + ui (—ms + uj)e fc 
(— tos + cj)e l * — v — u) y yd 



■(7) 



The system in this limit is interesting because it pro- 
vides a model to show the nonlinear effect when the 
nonlinearity appears only in the off-diagonal term. The 
self-trapping and the adiabatic evolution alike in the 
sense that both of them can conserve the population 
on a state. The difference is that the (traditional) self- 
trapping happens at a time-independent state, while the 
adiabatic evolution conserves the population on the time- 
dependent eigenstates of the Hamiltonian. We now show 
that the self-trapping can appear at a time-dependent 
state. By the adiabatic theorem, if the system starts in 
an eigenstate of the Hamiltonian at time t = 0, it will 
evolve to the corresponding eigenstate of the Hamilto- 
nian at later time as long as the Hamiltonian changes 
slowly enough. For the considered system without non- 
linearity (to = 0), the adiabatic condition [28| reads 



A = lu/2v < 1, 



(8) 



where A is the adiabatic condition that we will refer to in 
later discussion. From Eq. ©, we find that the change 

rate (£L| ) of the Hamiltonian affects the number and 

behavior of fixed points, now we study this effect by the 
phase space analysis. Define, 



z = Z= \d\ 2 - |c| 2 

e = e-^ = o d - 



where z and 9 satisfy 



i = — u)\l — z 2 sin( 



6 : 



■■ cos — mz cos 



f (1 - z 2 )sin26>, 
2 a t „, 



(9) 



(10) 



Then we can cast the dynamical system into a classical 
Hamiltonian 



H e (z 7 8) =vz— yz 2 cos 2 8 + -^ cos2#-wVl - z 2 cos( 

4 (11) 
By considering 9 — and i = 0, we obtain the following 
equation, 

z 4 -2^3_ (1 „^„^ )z2 + 2 - z _^l = . (12) 

TO TO TO TO TO 

To study the stability of the fixed points, we introduce 
the Jacobian matrix 29, 30] 



J = 



/ dz/dz 
\ 39 fdz 



dz/89 \ 
89/89 ) 



(13) 



It is well known that the eigenvalues of the Jacobian ma- 
trix characterize the type of fixed points and determine 
the stability of them. Two imaginary eigenvalues indicate 
a stable elliptic fixed point, while two real eigenvalues 
correspond to an unstable hyperbolic fixed point. So the 
stable elliptic fixed points can be singled out by calculat- 
ing eigenvalues of the corresponding Jacobian matrix. 
From Eq. (fT2"|) , we find that the number of the fixed 
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FIG. 3: The pole fixed point versus the nonlinearity m and cj, 
when v is fixed, oj is proportional to the adiabatic condition 
A. 



points depends on the ratio of the nonlinear parameter 

Mi to the adiabatic term When nonlinearity term 

is smaller than the adiabatic term, i.e., \M%\ < \E^ |, 

there are only two fixed points, while for strong nonlin- 
earity (| Ma | > | -Eh |), the other two fixed points appear 

in the phase space with proper coupling constant v, this 
was shown in Fig. [5] and in Fig. 0Ja). This is interesting, 

as E.\ in linear system can induce population transfer 

between states \E + ) and \EJ), while it might lead to bi- 
furcation in nonlinear systems. In what follows, we focus 
on the elliptic fixed point nearest to z = ±1. This fixed 
point is stable and describes approximately the instanta- 
neous eigenstate of the system, hence we call it pole fixed 
point. With strong nonlinearity to, the pole fixed point 
is very close to z = ±1, see Fig. [2] and Fig. [3] In this 
sense we claim that nonlinearity favors the fixed points 

2 = ±1. 

When 4 > 1, i.e., the adiabatic condition in linear 
case has been completely broken down, Eq. (fT2j) can be 
solved analytically up to the zeroth-order in v/m (we 
assume that to ^> v), there exist three solutions given by 



±\ 1- 



(14) 







Eq. (jT4"]) tells that we can manipulate the pole fixed point 
by adjusting the nonlinearity to to satisfy to 3> uj in this 
limit, this is numerically confirmed as shown in Fig. EJc). 



IV. ASYMMETRIC DOUBLE- WELL 
POTENTIAL 



N 




N 0- 



2 4 6 

t 



FIG. 4: (color online) (a) The distribution of the fixed points 
with different nonlinearity. The red-solid line depicts the 
pole fixed point and blue-dashed lines denote the other fixed 
points. Parameters chosen are oj = 1 and v = 1. (b) Time 
evolution of the population imbalance with eigenstates of 
the Hamiltonian as the initial state. Parameters chosen are 
oj = 1, v = 1. m — 0,2,10 correspond to blue-dash-dotted 
line, green-dashed line and red-solid line, respectively, (c) is 
the same as (b), but with different parameters, i.e., oj = 1, 
v — 0.01 and m — 0, 100, 400. The system is trapped ap- 
proximately in one of the instantaneous eigenstates even if 
adiabatic condition is broken, as the red-solid lines in (b) and 
(c) show. 



case, z and 9 satisfy 

UJV 



■VT 



-z\J 1 — z 2 sin( 



-z + 



mv 

UJV z 



a VT - z 2 



(1 - z 2 ) sin 29, 
cos 9 



mrv 1 — 2z 2 
+ — — cost 
a z VI - z 2 



-zcos 



(15) 



The effective Hamiltonian and the equation for z can be 
analytically expressed as 



m 

UJV 

a 
mv 2 
Act 2 



z 2 (r 2 -v' z cos^ 9) + 



mrv 



VT 



VT — z 2 cos ( 



cos 29, 



(16) 



In this section, we discuss a more general case that the 
BECs are trapped in an asymmetric double- well potential 
(r ^ 0). With the same definition for z and 9 in Eq. © , 
we map the system into a classical phase space. In this 



„,2 44, A/ 2 i 2 2 r> \ 2 

m a z + a [a + oj — to — 2ujr)z 

3/2 2 \3 / 3 22 22\ 

+2ma (r — v — ujr)z + 2ma(ujr — a r + a v )z 
+m 2 r 2 v 2 - a 2 (a 2 - uur) 2 = 0. (17) 
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FIG. 5: (color online) (a) Illustration of the biggest-valued 
fixed points (including stable and unstable ones) with differ- 
ent nonlinearity. (b) The biggest-valued stable elliptic fixed 
points (pole fixed point) with different nonlinearity. Here 
l x = \/l — z 2 cos (3, l y — \/l — z 2 sin/3 and l z = z, where z 
is the population imbalance. Parameters chosen are uj — 2, 
a — I and m — 0, 2, 4 corresponding to blue-dashed, green- 
dotted and red-solid lines. 




A 

FIG. 6: (color online) The pole fixed point versus /3(see (a)) 
and A (the adiabatic condition parameter, see (b)). f3 is in 
units of 7T. Parameters chosen are u) = 2, a = 1. m = 0, 2, 4 
are for the blue-dashed, green-dash-dotted and red-solid lines 
respectively. 



By the phase space analysis, we find that, different 
from the symmetric case, the fixed points are determined 
not only by the nonlinearity m but also by the ratio of 
the level bias r to the coupling constant v. Without 
nonlinearity, the adiabatic condition [28[ is 

A = ^- sin 8 < 1. (18) 
2a 



Now we investigate the dynamics of the system un- 
der different nonlinearity m and /3, which have the same 
notations as before, while u and a are fixed correspond- 
ing to a constant frequency and amplitude of the driving 
field. The adiabatic condition becomes A — sin/3 <C 1, 
when we set uj/2a = 1. By examining the location of 
the biggest-valued fixed points with different nonlinear- 
ity from adiabatic regime to diabatic regime, we find that 
the nonlinearity in BEC does play a positive role in the 
process of adiabatic evolution as shown in Fig. [5] and 
Fig. [5] When nonlinearity is weak, it pushes the pole 
fixed point to z = ±1 at all circumstances, as shown by 
the green-dotted line in Fig. [5jb). When the nonlinear- 
ity becomes strong, bifurcation appears, corresponding 
to the occurrence of new fixed points in the phase space 
and the system can evolve on its instantaneous energy 
eigenstate beyond adiabatic condition at some special /3, 
as shown by the red-solid line in Fig. [(^a). Before clos- 
ing this section, we would like to address that all fixed 
points can be easily calculated and analyzed. Aiming to 
study the self-trapping in the instantaneous eigenstates 
of the Hamiltonian, we here mainly focus on the fixed 
point nearest to z = ±1, discussions on the other fixed 
points can be carried out in the same way. 

Two remarks are now in order. (1) The linear and 
nonlinear systems discussed here are not independent, 
the linear system is exactly a limiting case of the nonlin- 
ear system with zero nonlinearity. (2) The fixed points 
are defined based on the basis spanned by the instan- 
taneous eigenstates, so the imbalance of the population 
represented by z denotes the population difference on the 
two instantaneous eigenstates, it ranges from —1 to 1 de- 
pending on the system parameters. Moreover, the popu- 
lation on the right and left wells can also be manipulated, 
this can be found via the definition of the instantaneous 
eigenstates. 



V. CONCLUSION 

In summary, the time-dependent self-trapping for 
Bose-Einstein condensates in a double- well potential has 
been introduced and studied in this paper. Both the 
atom-atom coupling and the quality which character- 
ize the change of the system play important roles in 
the self-trapping. Fixed points are calculated and dis- 
cussed. These results suggest that the nonlinearity can 
help tracking the instantaneous eigenstates of the time- 
dependent Hamiltonian, providing a way to manipulate 
quantum systems. 
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